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We present the results of an analysis of a 2+1 dimensional pure SU(A0 Yang-Mills theory formu- 
lated on a 2-dimensional spatial torus with non-trivial magnetic flux. We focus on investigating 
the dependence of the electric-flux spectrum, extracted from Polyakov loop correlators, with the 
spatial size /, the number of colours N, and the magnetic flux m. The size of the torus acts a pa- 
rameter that allows to control the onset of non-perturbative effects. In the small volume regime, 
where perturbation theory holds, we derive the one-loop self-energy correction to the single-gluon 
spectrum, for arbitrary N and m. We discuss the transition from small to large volumes that has 
been investigated by means of Monte-Carlo simulations. We argue that the energy of electric 
flux e, for the lowest gluon momentum, depends solely on e/N and the dimensionless variable 
x = AM, with A = g 2 N the 't Hooft coupling. The variable x can be interpreted as the dimension- 
less 't Hooft coupling for an effective box size given by Nl. This implies a version of reduction 
that allows to trade / by N without modifying the electric-flux energy. 
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1. Introduction 

In this paper we will explore the electric-flux spectrum of a 2+1 -dimensional SU(N) Yang- 
Mills theory defined on a 2-dimensional spatial torus endowed with a chromo-magnetic flux. The 
aim is to disentangle the dependence on the number of colours N and the size of the torus I, fol- 
lowing the idea of reduction introduced by Eguchi and Kawai [Jl], The size of the torus will 
act as a control parameter that permits to explore the onset of non-perturbative dynamics. This 
is so, because for small / the effective coupling constant becomes small and perturbation theory 
holds. As / grows, the finite size effects, including the magnetic flux, should become irrelevant. 
The way in which this takes place, and in particular the interplay with the large N limit, might 
shed some light into the processes relevant for non-perturbative dynamics. We will focus here for 
simplicity on the 2+1 dimensional case that shares many of the non-perturbative properties of the 
4-dimensional theory. There is an extensive literature on the subject of Yang-Mills 3-d fields and 
large ,/V dynamics, for recent lattice reviews we refer the reader to 

2. Set-up and Perturbative Analysis 

We will be considering SU(A0 Yang-Mills theories defined on a spatial torus of size I x I. In 
the basis of constant transition matrices T,-, the gauge field connection has to satisfy the periodicity 
condition: Aj(x + lj) = r ; A,(^;)rt , with the T, fulfilling: 

r 1 r 2 = ex P {/^}r 2 r 1 , (2.1) 

where m is the magnetic flux. In the case that m and ,/V are co-prime, this equation defines the 
matrices r, uniquely modulo global gauge transformations. For simplicity we will assume that ,/V 
is odd and co-prime with m in the following. The periodicity constraint on the gauge fields can be 
solved by introducing a basis of N xN matrices f (p c ) satisfying: 

r,-f(p c )rj = ^f(p c ) , (2.2) 

where p c = jff-, with integers defined modulo N. Thus, there are N 2 such matrices. An explicit 
solution to the equation is: 

t{p c ) = e ia{ p C) ry*" 2 rf ni , (2.3) 

where k is an integer satisfying mk = 1 mod N, and a(p c ) an arbitrary phase factor. In this basis 
we can Fourier decompose the gauge connection as: 

A i {x)=^Y f e ip *A i {P)t{P) . (2.4) 

P 

The total momentum P is quantized in units of 2n/{Nl) and can be decomposed as: P = p + p c , 
with p = 2nn/l, and n\ G 7L. The prime means that p c = is excluded from the sum, since A,(x) 
has to be traceless. This can be interpreted by saying that the total momentum P is composed of 
two pieces: a colour-momentum part p c and a spatial-momentum part p. If we neglect the prime 
in the sum, the range of values of P is just that of a theory defined on a box of size Nl x Nl. The 



2 



Volume dependence in 2+1 Yang-Mills theory 



Margarita Garcia Perez 



e = (N 9, 0) 



e = (N e, N e 









_ Lattice +' 
Continuum 








I 
+ 











~ Lattice 
Continuum 




Figure 1: We show, for electric flux e = (N9,0) (left), and e = (NO, NO) (right), the function G(8) that 
gives the one-loop correction to the energy of a one-gluon state, through Eq. (2.5). The blue line corresponds 
to Eq. (2.6), while the red points are derived using a lattice regularization in the calculation of the self-energy. 



formalism looks hence as if the theory had no colour, but was defined on a bigger spatial box with 
effective size Nl. 

We will now analyze in perturbation theory the energy spectrum of one-gluon states. The first 
remark is that, since zero momentum is excluded, the theory has a mass-gap even for small Nl. 
At tree level in perturbation theory, the mass of a state with colour momentum p L is Eo = \p c \ = 
2n\n\/ (Nl). Such states can be characterized in a gauge-invariant way in terms of the electric flux 
e, which is defined modulo N and related to n through: <?, = Sjitijk. The gauge invariant operator 
projecting onto a state of electric flux e is a Polyakov loop winding times around each of the 
cycles of the 2-torus. The non-perturbative spectrum of such states can hence be obtained from 
Polyakov loop correlators. 

In perturbation theory the energy of electric flux can be determined by computing the gluon 
self-energy on the twisted box (for examples in SU(2) see Refs. [Q, ||]). Details of the calculation 
at one-loop order will be presented elsewhere [||]. Here we simply provide the final expression that 
will be used to analyze the numerical results in section ||. Putting everything together, we obtain 
for the square of the energy of electric flux e with colour momentum n: 



X 2 V AM J XNl \N 

where X = g 2 N is the 't Hooft coupling, which is dimensionful in 2+1 dimensions. The first order 
correction to the tree level expression is given in terms of: 

with the Jacobi Theta function given by: 

03(z,ix) = exp{-*7Dz 2 +2ninz} ■ (2.7) 

The dependence of G on e/N is exhibited in Fig. [jl for electric fluxes of the form e = (N6,0) and 
e = (NO, NO). 
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One important remark at this stage is that, once n and e/N have been fixed, the energy de- 
pends exclusively on the variable x = XNl, which can be interpreted as the dimensionless 't Hooft 
coupling corresponding to an effective box of size Nl. Based on our numerical results, we will 
argue that this is the case even in the non-perturbative regime. This will imply thus a kind of re- 
duction in which N and / become interchangeable. A second remark in order is the non-analyticity 
of G(6) at 6 = 0(mod 1). This could give rise to a tachyonic instability in the energy of electric 
flux for N — > oo, which has been previously discussed in the context of non-commutative geometry 
in Ref. [Qj. There it was presented as an evidence of spontaneous symmetry breaking. We will 
come back to this issue when discussing the numerical results in section ^ 

Let us finally mention that in order to arrive at the one-loop perturbative expression we have 
used three different approaches: to compute the Euclidean gluon self-energy in the continuum and 
on the lattice, and to use the Hamiltonian formulation. The lattice calculation has made use of the 
results in Refs. [Q],[j8|]. In Fig. [j], the lattice results for G(6), extrapolated to the continuum limit, 



are compared with the formula given in Eq. ( |2.6[ ). There is perfect agreement between the two. 

3. Numerical Results 

In this section we will present the results of a lattice calculation of the electric-flux spectrum 
in a twisted box as a function of N, I and the magnetic flux. We start with a Lq x L 2 lattice with 
twisted boundary conditions (/ = La, with a the lattice spacing). It is possible to perform a change 
of variables of the link matrices that allows to work with periodic links at the price of introducing 
a twist dependent factor in the action: 

S w =Nb'£ 1 £ (tf-z* v (x)TrP M v(je)) , (3.1) 

where Zn V (x) = 1 except for the corner plaquettes in each (1,2) plane where it takes the value: 

f 27Tm -i 

Zij (x) = exp | i —— Stj | , (3.2) 

with m the magnetic flux. We have performed Monte Carlo simulations with this action employing 
a heat-bath algorithm based on a proposal by Fabricius and Haan [[K]]. Four sets of (N,L,Lq) values 
have been generated: (5,14,48), (5,22,72), (7,10,32), and (17,4,32), with values of the magnetic 
flux: m = 1, 2, for SU(5), m = 1, 2, 3, for SU(7), and m = 1, 3, 5, 8, for SU(17). 

The electric-flux spectrum has been extracted from spatially smeared Polyakov loop correla- 
tors. The spatial Polyakov loops, winding e times around the ?th direction of the torus, read: 

P i (t,x,e)=Tr{Y[Ui(t,x + si)Y . (3.3) 

L .v=0 J 

The twist gives rise to non-trivial boundary conditions for the loops: P{(t, x + Lj, e) = zfj Pi(t, x, e), 
that have to be considered when extracting the spectrum from the Polyakov loop correlators. In 
practice, we project over the different colour-momentum states by computing the averaged corre- 
lators: 

C l (t,e,n) = j ? ^ex V [-i—yj(P 1 (0,y%e)Pl(t,(y+y)%e)) , (3.4) 
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Figure 2: For different values of /, N, and e, we display the quantity 8E 2 (e)/X 2 = (E 2 — Eq)/X 2 vs 
x = AM, where E is the energy of electrix flux e with momentum \n\ = 1, and Eq = 2 asinh(sin(7r/ (NL))) 
is the lattice tree-level perturbative result. 



with e = £ 2 i nk. An analogous expression holds for Ci{t-, e, n). The energy is extracted from the 
exponential decay at large t of the correlator. 

A detailed account of the results will be presented in [Eh. Here we will focus on a few examples 
that illustrate the scaling of the electric-flux energies with x = XNl. We will focus on the minimal 
momentum states with \n\ = 1 which have electric flux \e\ =k. According to our perturbative 
calculation, for given e/N we expect energies to depend solely on x. Given that N is prime, identical 
ratios of e/N are not possible, nevertheless our data falls in three sets with similar values for 
the ratio: (N,e) = {(5, 1), (17,3)}; {(7,2), (17,5)}; {(5,2), (7,3)}. The dependence of the energy 
squared on x is presented in Fig. |2| The scaling holds in all the range of x analyzed which goes far 
beyond the perturbative regime. A discussion of the results will be presented below. 

4. Discussion 

Our numerical results indicate that the energy of electric flux, in units of the 't Hooft coupling, 
depends solely on the variable x = XNl and the value of e/N, at least for the minimal momentum 
\n\ = 1. This does indeed hold for the perturbative expression. Let us look now at the expectation 
in the large x regime, where strings of electric flux are expected to be formed and the energy 
should grow linearly with the box size /. A large amount of lattice results hint in the direction 
that flux tubes can be described in terms of an effective string picture based on the Nambu-Goto 
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x=mii x=m\ x=xn\ 

Figure 3: In the upper panels we display several fits of 8E 2 {e)/X 2 to Eq. (4.4). The symbol codes are 
the same as in Fig. 2. Different lines correspond to: (combined) a = G(e/N), j5 and y free; (1-loop) 
a = G(e/N), j3 = y = 0; (string) a — 0, j3 and /free. The lower panels display 8E 2 (e)x 2 /X 2 . 



action ph. In the presence of a constant background magnetic field B, the Nambu-Goto prediction 
for a closed-string winding e times around the torus is given by: 



A 2 \X 2 N J v v 3A 2 V V A/ 

where a is the fundamental string tension. The last term on the right hand side of this formula 
can be easily mapped onto the tree-level perturbative expression Eq by taking into account the 
relation between n and e, and identifying the magnetic field with B = 2nm/N. If the Nambu-Goto 
expression would hold, two different regimes would take place in terms of the scaling variable x: 



Low x, where: 



Large x, where: 



X 2 A 2 x \N 



Z^(£)_£^(£) ( o\e[\ 2 2 _na_ 
A 2 A 2 VA 2 A^/ 3A 2 ' ^ ; 

This would be fully consistent with our hypothesis on the x scaling, once e/N is fixed. A formula 
that captures the asymptotic behaviour in both the low and large x regimes is: 

8E 2 (e) E 2 {e) E 2 {e) Ana Kp , 2 
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Fitting this expression to our data gives the results displayed in Fig. [| The fits are restricted to the 
large x region with x > 20 (x > 30, for N =\1, e = 3). Two different fits are presented: (a) fixing 
a = G(e/N), denoted by 'combined', and (b) fixing a = 0, denoted by 'string'. In all cases, the 
'combined' fit has a slightly better % 2 per degree of freedom than the 'string' fit. In addition, we 
display the 1-loop expression. Although the combined formula fails to reproduce the results in the 
intermediate x region it describes qualitatively quite well the data. One remarkable observation is 
that for small x it improves the one-loop result in the correct direction. This is clearly observed 
in the low panels of Fig. || where 8E 2 x 2 /A 2 is plotted to enhance this effect. We are at present 
exploring alternative parametrizations to improve the fit in the intermediate x regime. 

One final remark refers to the possible existence of tachyonic instabilities. One could use 
Eq. (4.4) to see in which instances the energy remains non-tachyonic for all values of x. A more 
thorough analysis will be presented in Let us just mention here that our results indicate that it 
suffices to keep both m and k of order N to avoid the tachyonic behaviour. This coincides with the 
criteria introduced in ^ for twisted Eguchi-Kawai reduction to hold at large N. 



Acknowledgments 

M.G.P. and A.G-A acknowledge financial support from the MCINN grants FPA2009-08785 
and FPA2009-09017, the Comunidad Autonoma de Madrid under the program HEPHACOS S2009/ESP- 
1473, and the European Union under Grant Agreement number PITN-GA-2009-238353 (ITN 
STRONGnet). They participate in the Consolider-Ingenio 2010 CPAN (CSD2007-00042). M. 
O. is supported by the Japanese MEXT grant No 23540310. We acknowledge the use of the IFT 
clusters for our numerical results. 



References 

[1] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48, 1063 (1982). 

[2] A. Gonzalez-Arroyo and M. Okawa, JHEP 1007, 043 (2010); The string tension from smeared Wilson 
loops at large N, arXiv: 1206.0049 [hep-th]; PoS(Lattice 2012)221. 

[3] R. Narayanan and H. Neuberger, PoS LAT 2007 (2007) 020. M. Teper, Acta Phys. Polon. B 40 (2009) 
3249. B. Lucini and M. Panero, SU(N) gauge theories at large N, arXiv: 1210.4997 [hep-th]. 

[4] A. Gonzalez Arroyo and C. P. Korthals Altes, Nucl. Phys. B 311 (1988) 433. 

[5] D. Daniel, A. Gonzalez-Arroyo, C. P. Korthals Altes and B. Soderberg, Phys. Lett. B 221 (1989) 136; 
D. Daniel, A. Gonzalez-Arroyo and C. P. Korthals Altes, Phys. Lett. B 251, 559 (1990). 

[6] M. Garcia Perez, A. Gonzalez Arroyo and M. Okawa, Small volume Physics in the 2+1 dimensional 
Yang-Mills theory, to appear. 

[7] J. R. Snippe, Phys. Lett. B 389, 119 (1996); Nucl. Phys. B 498, 347 (1997). 

[8] M. Luscher and P. Weisz, [Erratum-ibid. 98 (1985) 433]; Phys. Lett. B 158 (1985) 250; Nucl. Phys. B 
266(1986) 309. 

[9] Z. Guralnik, R. C. Helling, K. Landsteiner and E. Lopez, JHEP 0205, 025 (2002). 
[10] K. Fabricius and O. Haan, Phys. Lett. B 143 (1984) 459. 



7 



